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ABSTRACT ■ 

An  approximate  boundary  for  a  finite  sample  sequential  decision 
process  is  presented  (without  derivation).  By  means  of  Computer  Simu- 
lation, properties  of  the  process  are  checked;  the  power  of  the  test  is 
thus  determined  for  selected  alternative  hypotheses. 
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I.   INTRODUCTION 

The  main  purpose  of  this  paper  is  to  study  the  properties  of  a 
mathematically  derived  sequential  sampling  procedure.   The  derivation 
involves  the  assumption  that  a  certain  parameter  N,  identified  as  a 
sample  size,  was  "large."  Our  purpose  was  to  make  computer  calculations 
that  would  reveal  how  the  sampling  procedure  might  behave  for  finite, 
or  even  rather  small,  N.  Other  calculations  would  also  indicate  the 
power  characteristics  of  the  test,  that  is,  the  probability  of  rejecting 
a  null  hypothesis.  The  null  hypothesis,  in  this  case,  was  that  a  partic- 
ular probability  distribution  describes  the  sample.  By  experimental 
sampling  methods  we  attempted  to  estimate  the  power  that  the  procedure 
possessed  to  detect  departures  of  various  sorts  from  the  null  hypothesis. 

The  stochastic  process  used  to  obtain  the  sequential  procedure  might 
also  be  used  for  other  purposes  in  operations  research.   For  example, 
the  sequential  procedure  might  be  used  to  describe  arrivals  in  transitory 
qu-.ueing  situations.  Some  of  these  situations  are  briefly  described  at 
the  end  of  this  paper. 


II.  A  SEQUENTIAL  SAMPLING  PROBLEM 

Suppose ,  for  example,  that  a  supplier  of  transistors  asserted  that 
each  of  his  transistors  had  an  in-service  lifetime  T,  where  T  was  a 
random  variable  distributed  exponentially; 

f  1  -  e"  Xol  9   t  *  o 
F(t)  =\ 

to        ,  t  <  0  . 

In  order  to  test  that  assertion  or  hypothesis,  a  sample  of  N  tran- 
sistors was  selected  and  placed  in  service  simultaneously  at  time  t  =  0. 
On  the  average,  the  number  of  transistors  that  failed  by  time  t  >0  was 
NF(t),  i.e.-  if  A.T(t)  denoted  the  number  of  transistors  that  failed  by 
time  t,  then 

E  [AN(t)3  =  NF(t) 

and 
VarJAN(t)]  =  NF(t)  [l-F(t)]  , 

since  the  transistors  were  presumed  to  fail,  or  die,  independently. 

If  one  waited  until  all  N  transistors  had  failed,  the  rate  could  be 
estimated  and  compared  to  -A  .  A  hypothesis  test  can  be  conducted  by 
waiting  until  r  out  of  N  transistors  failed  (r<N);  this  was  shown  by 
Epstein  and  Sobel  [1]  . 

The  present  approach  to  this  problem  was  to  establish  a  boundary, 
B^t),  with  the  property  that  if   >=  X  ,  i.e.  if  the  supplier  was  en- 
tirely honest,  then  Ajj(t)<  Bjj(t)  for  all  t,  with  a  known,  specifiable, 
probability  denoted  by  1-6.   For  instance,  suppose  6  =   0.0£.   Then,  if 
^  =  ^o>   ^he   Pr°bability  that  the  accumulated  numbers  of  failures  by 
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time  t  would  ever  exceed  the  boundary  is  five  per  cent.  Hence,  if  A.. 
remains  below  the  boundary,   A£  A  ,  with  ninety-five  per  cent  confidence 

The  detailed  derivation  of  the  boundary  is  given  by  Gaver,  Lehoczky, 
and  Perlas  in  an  unpublished  technical  report  that  is  in  preparation.  The 
result  is  quoted  below. 

The  Boundary 

If  F(t)  is  sufficiently  smooth  (possesses  a  continuous  density)  and 
if 

BN(t).=  NF(t)  +VN  f*F(t)  +p  Cl-F(t)]}  ,  (2.1) 

then  the  probability  that  Aft)  <  B  (t)  for  all  t  is  approximately 

-2*p 

1 -e     ;   the  approximation  improving  as  N  becomes  large  (N  — *  CD  ) .  Here, 

oi    and  p  are  arbitrary  positive  numbers . 

A  graph  of  a  possible  failure  arrival  pattern  is  a  step  function  with 
unit  increases  at  time  points  that  are  the  order  statistics  of  F(t): 
T/.J  \,  T(2)>    •••  where  T/^  \  is  the  smallest  of  the  N  failure  times,  T/?n 
is  the  next  smallest,  ...  T/,\  is  the  kth  smallest.  Then,  a  particular 
sample  may  be  represented  as  follows: 

A^(t)  =  0   for   0  i  t  <  T,  , 

«  1   for   I(1)*t<T(2) 

■  2  ?°r       T(2)  $t  <  T(3)        (2.2) 


-  k   for   T(k)<t  <T(k+1) 


-  N   for   I    <  t  , 


and  graphically,  we  have  the  following  picture. 


T(D  T(2)  T(3) 


Figure  1  . 
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In  Figure  1,  a  crossing  of  the  boundary  occurs  at  T  .  We  would  conclude 
that  the  null  hypothesis  is  not  satisfied  if  such  an  event  occurs.  Notice 
that  the  parameters  ot  and  "p  may  be  selected  at  will.  It  is  interesting 


to  arrange  that  c<  f>  =  k,   a  constant  such  that  e 


-2«f    =  -2k 


0.05, 


for  example,  and  then  to  attempt  to  select  ex  (  or  "^  )  so  that  the 
boundary  will  be  crossed  soon  if  a  non-null  distribution  of  specific 
sort  is  to  be  detected  [2j  .  This  problem  has  not  been  extensively  in- 
vestigated, but  the  expected  time  to  hit  the  boundary  has  been  tabulated 
for  certain  of  the  problems  studied. 


III.  USE  OF  THE  COMPUTER  SIMULATION 

Several  facts  about  the  boundary  %(t)  presented  in  equation  (2.1  ) 
are  worthy  of  notice . 

a.  The  prescribed  probability  of  crossing  the  boundary,  e    '  , 
is  not  exact  for  finite  N.  That  is,  use  of  this  probability  to  assess 
the  probability  of  crossing  is  an  approximation  which  is,  in  fact,  quite 
analogous  to  the  use  of  the  normal  distribution  as  an  approximation  to 
the  distribution  of  a  sum  of  independent  random  variables.  The  mathema- 
tical derivation  does  not  indicate  whether  the  actual  probability  of 
crossing  is,  perhaps,  consistently  lower  or  higher  than  the  actual  true 
probability.  More  complex  mathematical  methods  may  provide  useful  in- 
formation of  this  type,  but  a  straightforward  approach  is  that  of  simu- 
lation or  synthetic  sampling.  By  this  means,  it  can  be  seen  hew  increases 
in  N  improve  the  quality  of  our  approximation,  as  is  suggested  by  the 
mathematical  derivation. 

b.  The  relationship  of  the  true  crossing  probability  tc  the 
approximate  crossing  probability,  e"2°<p  for  given  N,  may  possibly  be 
influenced  by  the  choice  of  *  and  v>   .  Setting  «  B  =  k,  where  k  is  a 
constant  delivering  a  desired  probability  of  crossing  (for  %%,   k&1.k9), 
«  (and  hence  "B  )  may  be  adjusted  so  as  to  obtain  good  agreement  between 
the  theory  (e    =  e    r  )  and  the  true  probability  for  various  values 
of  N —  in  particular  for  small  N.  In  this  paper,  an  attempt  was  made 

to  study  this  problem  by  simulation. 

In  order  to  study  the  properties  of  the  boundary  %(t),  the  following 
simulation  was  carried  out. 


A.   SIMULATION  LOGIC 

1  .  The  null  hypothesis  distribution  was  chosen  and  called  FgCt). 
The  formula  used  to  calculate  the  boundary  values  was  as  follows: 

B^t)  =  NFQ(t)  +  >TF  0*F0(t)  +p  [1-F0(t)]}  .        (3.1) 

An  exponential  distribution  is  the  null  hypothesis,  as  in  the  illustra- 
tive problem;  hence,  for  the  initial  calculations,  the  equation 

%(t)  =  N(l-e_t)  +lfN  {?+   («-?)  (l-e_t)}        (3.2) 

was  used. 

2.  A  sampling  distribution  was  selected.   Initially,  it  was  the 
null  hypothesis  distribution,  so  that  the  accuracy  of  the  theoretical 
estimate  of  the  probability  of  crossing  the  boundary  could  be  checked. 
Later,  other  distributions  were  tried  in  order  to  determine  the  power 
for  departures  from  the  exponential  distribution. 

3.  A  sample  of  size  N  was  taken  and  arranged  in  order: 

t(1)<t(2)<t(3)...<t(k)<...t(.0.  (3.3) 

This  ordered  sample  was  composed  of  the  actual  numbers  obtained  as 
realizations  of  the  random  variables. 

To)<T(2)<---T(k)<  ■"*<■>•  °-h) 

The  details  on  how  the  random  samples  were  conducted  will  be  presented 
shortly. 

k'     At  each  failure  time,  say  t/jx,  starting  with  k  =  1,  then  2,  3, 
...,  AN(t^k))  =  k  was  compared  with  %(t/k)):  if  k<%(t(k))j  for  k  = 
1,  2,  ...,  the  process  continued  to  k  +  1  and  was  checked  again.  If  at 
k*,  the  event  k*  £  ^^fk-'V  occurred  ?or   the  first  time/  a  boundary 


crossing  was  recorded  at  the  time  t  =  t^j  the  test  was  stopped  and 
another  simulation  begun. 

5.  Step  four,  above,  was  independently  repeated  NR  times  (NR  means 
the  number  of  replications).  The  probability  of  crossing  f^  was  esti- 
mated by  the  fraction 

number  of  crossings.  ,   ^s 

fN  ~         NR  UO; 

Clearly,  NRfN  was  a  realization  of  a  binomially  distributed  random 
variable.  An  approximate  idea  of  the  variability  of  the  estimate  of 
the  probability  of  crossing  is  given  by 

62    B_i L.  •  (3.6) 

NR       NR 

An  approximate  confidence  limit  for  the  true  probability  P  of  crossing 

2r  2 

came  from  assuming  that  f^  is  approximately  normal  with  variance     f.TP* 

B.   SIMULATION  PROGRAM 

The  simulation  procedure  just  described  was  written  in  FORTRAN  com- 
puter programming  language  and  run  on  the  IBM  3b0,  Model  6?  Computer  at 
the  Naval  Postgraduate  School.  The  program  was  written  to  facilitate 
changes  of  distributions,  boundary  parameters  ( <x  and  f>  ),  and  numbers 
of  replications  (NR).  It  was  also  constructed  so  that  the  sampling 
process  could  be  stopped  after  a  fixed  time  period,  i.e.-  before  all  of 
the  order  statistics  had  been  recorded. 

The  pseudo  random  number  generator  employed  to  obtain  realizations 
of  uniform  random  numbers  was  a  library  subroutine  called  RANDOM  f3  and  1$  ; 
it  is  based  on  a  Lehmer  multiplicative  scheme.   The  uniform  realizations 
were  then  transformed  into  other  random  numbers  by  a  probability  integral 
transformation . 
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A  program  listing  and  sample  output  has  been  included  at  the  end 
of  this  paper.  In  order  to  generate  the  sequence  of  order  statistics 
(interpreted  as  failure  or  arrival  times),  the  random  sample  of  size  N 
was  first  created.  Then,  a  SORT  routine  was  employed  to  put  the  times 
in  order.  A  more  efficient  procedure  would  have  been  to  generate  the 
ordered  increments  sequentially;  for  example,  the  time  between  T/,  \  and 

T/,  i  \  in  exponential  sampling  is  also  exponential  with  a  mean  propor- 

1 
tional  to  vr-jr-  One  program  has  been  written  in  this  vein,  but  production 

runs  for  the  present  investigation  have  largely  been  carried  out  by  use 

of  the  SORT  routine  C£J. 
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IV.   SIMULATION  EXPERIMENTS 

Several  simulation  experiments  have  been  carried  out  to  develop  an 
understanding  of  the  boundary  crossing  process. 

A.  EXPERIMENT  1 

This  experiment  explores  the  crossing  probability  under  the  null 
hypothesis  (  ^  =  ^0)  and  under  selected  alternatives,  for  various  sample 
sizes  (N)j  and  for  different  boundary  parameters  o(  and  ^  ,  but  with 
c<   =    f  .  The  computed  probability  of  a  crossing  fl.  was  calculated  by 
the  fraction  (3.5)  mentioned  previously.   The  average  number  of  arrivals 
up  to  a  crossing  (A~(t))  was  calculated  by  summing  the  number  of  arrivals 
(Ajr(t))  that  had  occurred  up  to  the  time  of  crossing,  over  the  two  thou- 
sand replications  and  dividing  by  the  number  of  crossings  (NC)  recorded. 
The  average  crossing  time  (C(t))  was  calculated  by  adding  the  arrival 
times,  which  generated  crossings,  then  dividing  that  sum  by  the  number 
of  crossings.  The  results  appear  in  Table  I  (theoretical  crossing  prob- 
ability of  five  per  cent),  and  Table  II  (theoretical  crossing  probability 
of  ten  per  cent).   These  Tables  are  based  on  two  thousand  replications 
(NR  =  2000). 

When  the  boundary  was  not  crossed  for  a  sample  size  N,  all  values  for 
that  N  appear  as  0 .  In  such  a  situation,  the  zero  (0)  is  intended  as  a 
meaningless  place  holder  rather  than  the  unique  number. 
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TABLE  I. 


5%  Probability  of 

a  Crossing 

«=  ■£  =  1  .22388 

EXPON   (    A  - 

.5) 

N 

NC 

A(t) 

c(t) 

% 

k 

3 

3-33 

.133 

0.15 

6 

2 

5-00 

.300 

0.10 

8 

7 

6.1*2 

•  U28 

0.35 

10 

1 

8.00 

.523 

0.05 

25 

0 

0.00 

.000 

0.00 

50 

0 

0.00 

.000 

0.00 

75 

0 

0.00 
EXPON  (    X  = 

1) 

.000 

0.00 

N 

NC 

A(t) 

c(t) 

fN 

h 

U6 

3-63 

.287 

2.30 

6 

69 

5-07 

.391 

3-1*5 

8 

66 

6.30 

•  U02 

3.30 

10 

70 

7-37 

.1*16 

3.50 

25 

73 

15.19 

.1*61 

3-65 

50 

•      85 

27.5b 

•50U 

I1.25 

75 

10lt 

U2-U3 

.586 

5.20 

EXPON   (    A-  1  .25) 


N 

NC 

A(t) 

c(t) 

% 

h 

100 

3.66 

.298 

5-00 

6 

11*1 

5.11 

.381 

7.05 

8 

165 

6.33 

.mil 

8.25 

10 

197 

7-1*8 

.1*28 

9.8U 

25 

372 

15.95 

.513 

18.59 

50 

61*1* 

28.61* 

•  51i2 

32.19 

75 

839 

1*1.09 

.560 

1*1.95 
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TABLE  I.      (continued) 
EXPON   (    X  =  1 .$) 


N 

NC 

A(t) 

c(t) 

Jl 

k 

17li 

3.68 

.303 

8.70 

6 

253 

5-03 

.362 

12.61* 

8 

331 

6.26 

•  395 

16.5U 

10 

iioU 

748 

426 

20.19 

25 

896 

15-67 

•  1*98 

Mi.  79 

50 

461 

25- 5U 

470 

73.  ok 

75 

1770 

35-59 

EXPON   (    A  =   1  .75) 

432 

88  49 

N 

NC 

A(t) 

c(t) 

fa 

k 

262 

3-70 

.301 

13.10 

6 

389 

1|.97 

.31*9 

1945 

8 

525 

6.19 

.381 

26.2U 

10 

672 

743 

418 

33.59 

25 

1391 

4.80 

437 

69. 5U 

50 

1893 

23.31 

.362 

9k -6h 

75 

1985 

28.30 

.27U 

99. 2U 

EXPOIJ   (    A  =  2 ) 


N 

NC 

A(t) 

c(t) 

i 

k 

360 

3.69 

.296 

1845 

6 

552 

ii-98 

.31*6 

27-59 

8 

7U1 

6.15 

.372 

37. oU 

10 

929 

7.31 

•  397 

1*6.14 

25 

1736 

13.95 

.382 

86.79 

50 

1985 

19.86 

.256 

99.2k 

75 

2000 

23-65 

.190 

100.00 
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TABLE  II. 

10%   Probability  of  a  Crossing 
*=  f   =  1  .07298 

EXPON  (  >>  =  .5) 


N      ■ 

NC 

A(t) 

c(t) 

i 

k 

15 

3-33 

.300 

0.75 

6 

9 

5.00 

-Il67 

o.ii5 

8 

11 

5.72 

•3U1 

0.55 

10 

5 

6.80 

•Uoo 

0.25 

25 

0 

0.00 

.000 

0.00 

50 

0 

0.00 

.000 

0.00 

75 

0 

0.00 
EXPON   (    A  =   1  ) 

.000 

0.00 

N 

NC 

A(t) 

c(t) 

ffi 

k 

127 

3.1+9 

.388 

6.35 

6 

151 

U-90 

.U28 

7.5U 

8 

I51i 

5.92 

.1x21 

7.70 

10 

156 

7.03 

•  iiU3 

7-79 

25 

160 

1U-86 

.ii95 

8.00 

50 

177 

26.7U 

.522 

8.85 

75 

190 

39.06 

.551 

9.50 

EXPON  (  ^  -  1  .25) 


N 

NC 

.     A(t) 

c(t) 

fk 

k 

2li1 

3.51 

.31+5 

12.05 

6 

289 

ii.72 

.381 

u.ii5 

8 

33U 

5.98 

.1*25  • 

16.70 

10 

371; 

7.10 

.1|2|6 

18.70 

25 

627 

15.U 

.517 

31.3a 

50 

9k$ 

26.86 

.527 

a7.2a 

75 

1168 

38.66 

•  51i3 

58.39 
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TABLE 

II.      (c 

ontinued ) 

EXPON   (    * 

-i-5). 

N 

NC 

A(t) 

c(t) 

fN 

h 

360 

3.51; 

.5h$ 

17.99 

6 

14.60 

It. 72 

.376 

22.79 

8 

571 

5.85 

.I405 

28.51* 

10 

671; 

7.01; 

-1+39 

33.69 

25 

1197 

1U.U5 

.hT2 

59.81; 

50 

1682 

21*.  00 

.1*21* 

81;. 09 

75 

1907 

31.37 

.372 

95-31; 

EXPON   (    >=  1  .75) 


N 

NC 

A(t) 

c(t) 

5l 

k 

1*98 

3.53 

.337 

21;. 89 

6 

660 

I*.  71* 

.381 

32.99 

8 

823 

5.81* 

.399 

1*1.11* 

10 

972 

6.88 

.1*10 

U8.59 

25 

1657 

13.59 

.1*11; 

82.81* 

50 

1957 

20.36 

.307 

97.81* 

75 

1998 

2l*.l*6 
EXPON    (    >  =  2 ) 

.230 

99.89 

N 

NC 

A(t) 

c(t) 

i 

h 

617 

3.51 

.326 

30.81* 

6 

868 

U-71 

.370 

1*3.39 

8 

1083 

5.77 

.387 

51;.  11* 

10 

121*7 

6.71; 

.386 

62.31* 

25 

1878 

12.1*0 

.337 

93.89 

50 

1997 

17.22 

.211* 

99.81* 

75 

2000 

20.51 

.160 

100.00 
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B.  DISCUSSION 

The  actual  probability  of  crossing  appeared  to  be  consistently  smaller 
than  the  theoretical  probability  under  the  null  hypothesis  (A  =  1  ).  For 
N  of  fifty,  and  over,  the  agreement  appeared  to  be  reasonably  good.  More- 
over, the  probability  of  rejecting  the  null  hypothesis  (power)  was  quite 
high  when  h  =  1  ,J>  and  N  was  at  least  twenty-five.  This  would  not  be  an 
unreasonable  sample  size  for  the  transistor  testing  illustration,  but 
would  undoubtedly  by  large  if  the  procedure  were  applied  to  check  entire 
sys  terns . 

Notice  that  an  empirically  generated  correction  could  be  derived  to 
adjust  the  size  of  the  null  hypothesis  probability.  Consider,  for  ex- 
ample, the  N  =  h}      A  =  1 .25  element  of  the  nominal  five  per  cent  table: 
the  crossing  per  cent  (estimated  true  probability  of  crossing)  is  nearly 
five  per  cent.  Under  the  null  hypothesis,  the  simulation  crossing  prob- 
ability was  computed  at  2.3  per  cent.  Now,  if  the  null  hypothesis  pre- 
vailed, and  if  each  failure  time  were  multiplied  by  (1 .25)   =0.8,  then 
five  percent  boundary  crossings  would  be  realized.  On  the  other  hand, 
if  A  >  1 ,  and  each  failure  time  was  multiplied  by  0.8,  but  still  applied 
against  the  theoretical  derived  boundary  for  A  =' 1 ,  the  power  of  the 
procedure  would  be  further  enhanced.   This  empirical  adjustment  has  been 
left  for  more  extensive  development,  e.g.  —  for  different  sample  sizes — 
in  the  future.  It  seems  clear  that  an  alternate  approach  involves  an 
adjustment  of  the  boundary,  itself. 

C.  EXPERIMENT  2. 

Using  the  same  procedure  as  in  Experiment  1 ,  we  explore  the  effect 
of  <*  j^   £  >   but  with  «*  f>    =   constant.  We  tried  the  following  four  com- 
binations:  (1  )  <*  =  1.5  P,   (2)  f   =  1.£*  ,   (3)  <*  =  2*p  ,  (h)   ?  =  2*  . 
The  results  are  summarized  in  the  next  tables. 
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TABLE  III 


5%  Probability  of  a  Crossing 
c*=  1.U9893    f  =   -99929 


EXPON   (    *  =.5) 

N 

NC 

A(t) 

c(t) 

f 

_N 

k 

2 

ll.OO 

.156 

.10 

6 

1 

6.00 

.1*76 

.05 

8 

3 

7.33 

•561* 

.15 

10 

0 

0.00 

.000 

.00 

25 

0 

0.00 

.000 

.00 

50 

0 

0.00 

.000 

.00 

75 

0 

0.00 

.000 

.00 

EXPON   (   >v  =  1  ) 


N 

NC 

A(t) 

c(t) 

fk 

li 

17 

ll.OO 

.297 

.85 

6 

36 

5.61 

.1|61 

1 .80 

8 

Mi 

6.95 

.1*92 

2.20 

10 

h9 

8.20 

.50a 

2.1*5 

25 

h9 

17.28 

.589 

2.1*5 

50 

70 

32.91 

.697 

3.50 

75 

91 

U7.U6 

.706 

U.55 

EXPON   (    X=  1  .25) 


N 

NC 

•  A(t) 

c(t) 

b. 

1* 

36 

ll.OO 

•  309 

1  .80 

6 

97 

$.$9 

.1*57 

1*.85 

8 

118 

7.00 

.1*93 

5.90 

10 

11*7 

8.31* 

.529 

7.35 

25 

31*0 

18.20 

.668 

16.99 

50 

636 

32.99 

.697 

31.79 

75 

871 

1*6.51* 

.692 

1*3.51* 
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TABLE 

III     (continued) 

EXPON   (    *  = 

i'.5) 

N 

NC 

A(t) 

c(t) 

fa 

h 

76 

U.oo 

.316 

3.80 

6 

161 

$•$9 

.1428 

8.0U 

8 

232 

7.03 

.h9$ 

11  .60 

10 

321* 

8.36 

.531 

16.19 

25 

857 

17.1*7 

.603 

1*2.81* 

50 

11*72 

29.77 

.562 

73.59 

75 

1799 

U0.3ii 

•  521 

89.91+ 

EXPON   (    >v  =  1  .75) 


N 

NC 

A(t) 

c(t) 

!n 

1* 

126 

li.00 

.31 U 

6.30 

6 

251* 

5.55 

.1*15 

12.70 

8 

U01 

6.91 

•  U70 

20. Ol| 

10 

560 

8.27 

.511 

27.99 

25 

1381 

16.67 

.538 

69.0k 

50 

1897 

26.  U2 

.1*36 

9l*.81* 

75 

1990 

32.57 
EXPON   (    A  -  2) 

•  33U 

99.  h9 

N 

NC 

A(t) 

c(t) 

i 

h 

189 

U.oo 

.310 

9.1U1 

6 

U02 

5.5o 

.1*09 

20.09 

8 

$9h 

6.85 

.k$k 

29.70 

10 

787 

8.11; 

.1*81 

39. 3U 

25 

1733 

15.70 

.1*62 

86.6U 

50 

1987 

22.91 

.312 

99.31* 

75 

2000 

27.50 

.229 

100.00 
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TABLE  IV 


$%   Probability  of  a  Crossing 
*  =  .99929    f  -  1 .U9893 


EXPON  (  >\  =   .5) 


N 

NC 

I(t) 

c(t) 

i 

k 

10 

3.10 

.203 

0.50 

6 

k 

1|.50 

.273 

0.20 

8 

10 

5.70 

.331 

0.50 

10 

3 

5.66 

.21*0 

0.15 

25 

0 

0.00 

.000 

0.00 

50 

0 

0.00 

.000 

0.00 

75 

0 

0.00 
EXPON   (   >  =  1 ) 

.000 

0.00 

N 

NC 

A(t) 

c(t) 

fa 

ii 

80 

3.30 

.2l|1 

ii.00 

6 

99 

u-Jua 

.301 

ll.95 

8 

90 

5.10* 

.291 

U.50 

10 

101 

6.27 

.30)4 

5.05 

25 

91 

12.68 

.331 

ii.^ 

50 

99 

23.08 

.371; 

i4.95 

75 

117 

3U.13 

-lj.11 

5.85 

EXPON  (  A  =  1  .25) 


N 

NC 

A(t) 

c(t) 

i 

II 

162 

3.38 

.263 

8.10 

6 

198 

1;.52 

.301 

9.90 

8 

21  It 

5. 61 

.323 

10.70 

10 

2U2 

6.57 

.331* 

12.10 

25 

392 

13.78 

.392 

19.59 

50 

628 

2U.36 

.1*13 

31.39 

75 

797 

35.16 

.1(31* 

39.8U 
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TABLE 

IV     (continued) 

EXPON   (    >  = 

Lfl 

N 

NC 

A(t) 

c(t) 

i 

h 

2l*8 

3.1*2 

.268 

12.39 

6 

33U 

h.hl 

.291 

16.70 

8 

1*00 

5.59 

.319 

19.99 

10 

htt 

6.58 

.333 

22.71* 

25 

889 

13.77 

.396 

W4.U5 

50 

11*06 

23.  1*1* 

.386 

70.29 

75 

1720 

31-21* 

.357 

85.99 

EXPON   (  >>  =  1  .75) 


N 

NC 

A(t) 

c(t) 

& 

i* 

31*6 

3.1*3 

.261; 

17.29 

6 

1*77 

1*.1*7 

.287 

23.81* 

8 

588 

$•& 

.313 

29.39 

10 

729 

6.60 

.335 

36.1*5 

25 

1368 

13.16 

.361 

68.39 

50 

1867 

20.58 

.303 

93.314 

75 

1979 

21*. 88 

.235 

98.95 

EXPON   (    X  =  2  ) 


N 

NC 

A(t) 

c(t) 

5l 

1* 

1+66 

3.1*5 

.268 

23.29 

6 

m 

U.51 

.293 

32.19 

8 

806 

5.53 

.307 

1*0.29 

10 

977 

6.53 

.323 

1*8.81* 

25 

1711* 

12.37 

.316 

85.70 

50 

1980 

17.39 

.216 

98.99 

75 

not  run 
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TABLE  V 


■\ 


5$  Probability  of  a  Crossing 
<X=  1.73082     p=   .8651*1 


EXPON  (  \=   .5) 


N 

NC 

A(t) 

c(t) 

£n 

k 

2 

I*.00 

.156 

0.10 

6 

1 

6.00 

.1*76 

0.05 

8 

1 

8.00 

.61$ 

0.05 

10 

0 

0.00 

.000 

0.00 

25 

0 

0.00 

.000 

0.00 

50 

0 

0.00 

.000 

0.00 

75 

0 

0.00 

.000 

0.00 

EXPON    (    h  =   1  ) 


N 

NC 

A(t) 

c(t) 

5l 

ll 

k 

it. 00 

.169 

0.20 

6 

17 

5.88 

.kh$ 

0.85 

8 

28 

7.50 

.$6h 

1  .JiO 

10 

3h 

8.76 

.567 

1  .70 

25 

39 

18.92 

.703 

1 .95 

50 

59 

36.86 

.858 

2.95 

75 

85 

52.07 

.81i5 

a. 25 

EXPON   (    *  =  1  .25) 


N 

NC 

A(t) 

c(t) 

fk 

l* 

11 

U.oo 

.207 

0.55 

6 

1*8 

5.87 

.1*60 

2.140 

8 

70 

7.1*0 

.538 

3.50 

10 

103 

8.81* 

.583 

5.15 

25 

308 

19.81* 

.805 

15.1*0 

^0 

612 

35.67 

.808 

30.59 

75 

851 

50.50 

.803 

1*2.51* 
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TABLE 

V     (continued) 

EXPON   (    )v  = 

i.5). 

N 

NC 

A(t) 

c(t) 

& 

1* 

17 

U.oo 

.198 

0.85 

6 

9k 

5.86 

•UU1 

l*-70 

8 

161 

7-iiO 

.531* 

8. Oil 

10 

253 

8.89 

^9$ 

12.6U 

25 

80i| 

18.82 

.691 

U0.19 

50 

11*63 

32.1*2 

.651 

73-11* 

75 

1802 

U3-71 

.589 

90.09 

EXPON   (  >  =  1  .75) 


N 

NC 

A(t) 

c(t) 

% 

h 

28 

U.oo 

.201* 

1  .1*0 

6 

169 

5.81 

.1*27 

8.U5 

8 

286 

7.33 

.507 

11*.29 

10 

1*1*7 

8.83 

•  571 

22.31* 

25 

1327 

18.03 

.620 

66.31* 

50 

1895 

28.73 

.1*92 

91*. 71* 

75 

1987 

35-92 

EXPON    (    X  =   2) 

.381 

99-3li 

N 

NC 

A(t) 

c(t) 

% 

h 

5U 

li.00 

.216 

2.70 

6 

239 

5.81 

.1*18 

11.95 

8 

1*1*7 

7.29 

.1*97 

22.31* 

10 

666 

8.70 

.538 

33.29 

25 

1721* 

17.12 

.535 

86.20 

50 

1985 

25-28 

.358 

99-21* 

75 

2000 

30.71* 

.265 

100.00 

23 


TABLE  VI 

$%   Probability  of  a  Crossing 
«=  .865U1  ?  =   1  .73082 


EXPON  (  *  =   .5) 


N 

NC 

A(t) 

c(t) 

% 

h 

111 

2.78 

.151 

0.70 

6 

6 

1*.16 

.258 

O.30 

8 

11 

$.$h 

.31  k 

0.55 

10 

5 

l*.80 

.176 

0.25 

25 

0 

0.00 

.000 

0.00 

50 

1 

13.00 

.130 

0.05 

75 

1 

18.00 
EXPON    (    X  =  1  ) 

.130 

0.05 

N 

NC 

A(t) 

c(t) 

% 

h 

97 

3.02 

.203 

U.85 

6 

117 

3.95 

.230 

5.85 

8 

106 

5.00 

.251* 

5.30 

10 

115 

5.60 

.239 

5-75 

25 

105 

11  .1*8 

.281 

5.25 

50 

119 

20.31 

-30l» 

5.95 

75 

122 

30.05 

.336 

6.10 

EXPON  (  >  =  1  .25) 


• 

N 

NC 

■  A(t) 

c(t) 

fJL 

il 

192 

3.10 

.225 

9.60 

6 

228 

U-13 

.253 

11  .1*0 

8 

21*5 

5.10 

.271* 

12.21* 

10 

265 

5.98 

.279 

13.25 

25 

1*00 

12.1*2 

.328 

19.99 

50 

600 

21.1*2 

.337 

29.99 

75 

777 

31  .81 

.373 

38.81* 

21* 


TABLE  VI     (continued) 
EXPON   (    *=Kj) 


N 

NC 

A(t) 

c(t) 

i 

1* 

299 

3.15 

.232 

ii*.  95 

6 

361 

l*.ii* 

.253 

18.01; 

8 

1*22 

5.08 

.269 

21.09 

10 

U81 

6.08 

.287 

2U.0U 

25 

859 

12.28 

•  323 

12.95 

50 

1363 

21  .33 

.331 

68.1  U 

75 

1679 

28.53 

.312 

83. 9U 

EXPON   (   >  =  1  .75) 


N 

NC 

A(t) 

c(t) 

3l 

h 

399 

3.15 

.226 

19.91* 

6 

501* 

U.l  7 

.251* 

25.20 

8 

613 

5.07 

.263 

30.6U 

10 

728 

6.n3 

.281 

36.39 

25 

1320 

11  .96 

.309 

65.99 

50 

18)43 

19.00 

.271 

92.11* 

75 

1972 

22.88 

.212 

98.59 

(  >i   ---   2) 


N 

NC 

A(t) 

c(t) 

5l 

1* 

526 

3.15 

.229 

26.29 

6 

682 

1|.15 

.250 

3U.09 

8 

837 

5.06 

.262 

1*1.81* 

10 

991 

6.01 

.280 

h9.$h 

25 

1703 

11.36 

.278 

85.11* 

SO 

1976 

15.98 

.191* 

98.79 

75 

not  run 

25 


D.  DISCUSSION 

The  tendency  of  the  theory  to  understate  the  crossing  probability, 
noted  <*.=  £  ,  is  even  more  pronounced  when  c<=  1  ,$f  .  However,  when 

"p  =  1  .5°^  (TABLE  IV)  the  understatement  error  nearly  vanishes,  and  the 
test  is  very  close  to  being  exact  (i.e.-  crossings  occur  with  probability 
about  five  percent)  for  all  values  of  N  checked.  According  to  TABLE  VI, 
the  probability  of  rejection  is  slightly  overstated  when  p  =  2°<v  ,  for 
all  but  M  =  k*     To  within  limits  of  sampling  accuracy,  it  appears  that 
"jj  equalling  1  .65°(  will  give  a  very  nearly  correct  answer  at  S   =  0.0S>. 
Further  sampling  will  pin  this  down  more  precisely.  A  theoretical  ex- 
planation would  be  most  interesting,  for  the  asymptotic  (large  N)  mathe- 
matics detects  no  difference  in  situations  for  which   <*fi  =  constant. 
The  next  set  of  experiments  was  conducted  in  order  to  investigate 
our  procedure's  ability  to  detect  departures  from  the  null  (exponential, 

A  =  1 )  hypothesis. 

E.  EXPERIMENT  3. 

Explorat:    ?   the  crossing  probability      :   boundary  reflected 
the  null  hypothesis  (exponential  with  ^  =  )*  q  =  1  ,  °^=pj  nominal 
crossing  probability  five  percent),  but  samples  were  drawn  from  a  special 
long-tailed  distribution  with  a  Pareto-like  appearance : 

cDC 

FP(x)  =  T^r'xJ0-  &•'> 

In  order  to  give  a  fair  comparison,  the  parameter  a  was  chosen  so  as  to 

match  the  medians  of  F  and  the  null  distribution;  thus  a  =  (log»2)   . 

P  e 

A  rationale  for  the  choice  of  (U.1  )  is  that  such  a  distribution 

appears  very  similar  to  the  exponential,  but  has  a  decreasing  failure 

rate  and  a  bias  toward  early  failures.  The  results  are  recorded  in 

TABLE  VII. 
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TABLE  VII 

5%   Probability  of  a  Crossing 
o<=  p  =  1.22388 


PARE  TO 


N 

NC 

A(t) 

c(t) 

i 

h 

79 

3.50 

.253 

3.95 

6 

91 

U.82 

.319 

U.85 

8 

91 

5. 61 

.2714 

U-55 

10 

1QU 

6.70 

.309 

5.20 

25 

133 

12.30 

.277 

6.6U 

50 

209 

20.00 

.257 

10.  \6 

75 

218 

26.65 

.2U2 

10.90 

F.  DISCUSSION. 

Our  sampling  experiment  indicates  that  the  sequential  procedure  under 
study  has  very  little  power  against  the  "Pareto  tail"  alternative  (1|.1  ) 
under  the  particular  choice  of  parameters ,  at  least  for  small  N.  Actually, 
the  probability  of  rejection  is  somewhat  larger  than  for  the  correspond- 
ing exponential.   It  seems  likely  that  the  relative  prevalence  of  small 
observations  accounts  for  the  increased  crossing  probability,  and  also 
for  the  relatively  short  expected  time  to  cross ,  given  that  a  crossing 
occurred. 

G.  EXPERIMENT  k- 

The  procedure  used  in  Experiment  3  was  repeated,  except  that  samples 
were  drawn  from  a  shorter-tailed  (Weibull,  with  increasing  hazard  rate) 
distribution: 

Fw(x)  =  1-e_ax2  .  U.2) 

Again,  the  parameter  a  was  fitted  by  matching  medians  of  the  exponential 
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-1 

to  (li.2);  by  coincidence,  a  =  (loge2)   once  again —  an  accident  of 
parameterization.  The  results  are  recorded  in  TABLE  VIII. 


TABLE  VIII 


$%   Probability  of  a  Crossing 
*=  p  =  1  .22388 


N  NC 


h 

13 

6 

38 

8 

58 

10 

91 

25 

805 

50 

1891; 

75 

1996 

A(t) 

c(t) 

i 

li.00 

.ll2S 

0.65 

5.89 

.610 

1  .90 

7.7.2 

.702 

2.90 

9.1|6 

.778 

h.55 

22.78 

1  .078 

U0.21* 

U1.88 

1  .091 

9ii.  69 

58.30 

1 .009 

99.79 

H.  DISCUSSION. 

The  sampling  experiment  suggests  that  the  power  of  our  sequential 
procedure  is  quite  unsatisfactory  against  the  particular  Weibull  alterna- 
tive examined  for  small  N,  but  the  power  rapidly  increases  with  N.   Scru- 
tinizing of  the  density  of  (I4..2)  indicates  that  there  are  relatively  few 
"early  failuros",  while  a  concentration  of  "wearouts"  probably  accounts 
for  the  rather  long  conditional  expectation  of  the  time  to  hit  the 
boundary  (reject  the  null  hypothesis). 
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